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Abstract 

By applying a standard solution-generating transformation to an arbi- 
trary vacuum Bianchi type II solution, one generates a new solution with 
spikes commonly observed in numerical simulations. It is conjectured that 
the spike solutions are part of the generalized Mixmaster attractor. 

1 Introduction 

Berger and Moncrief [IJ studied Gowdy spacetimes and found small-scale spa- 
tial structures develop on approach to the initial singularity. Since then many 
efforts have been spent trying to understand these spiky structures in Gowdy 
spacetimes and in more general G2 spacetimes through numerical simulations 
and analytical approximations [H [31 [H O HI [7] ■ 

My motivation in studying spikes is to understand its role on approach to 
generic singularities. Lifshitz, Khalatnikov and Belinskii [5l[Hl[TD] were the first 
to provide heuristic arguments that the approach to generic spacelike singular- 
ities are vacuum dominated, local, and oscillatory (known as the BKL conjec- 
ture). Further evidence came from the study of Bianchi type IX cosmologies by 
Misner [TTJ[Ill[I3]j who coined the term "Mixmaster" to describe the oscillatory 
behaviour. Uggla et al [IT provided a detailed description of the local attractor 
for generic singularities (called the generalized Mixmaster attractor). See [15j 
for a more complete introduction and the latest work on the attractor (called 
the billiard attractor). The local part of the BKL conjecture is increasingly 
under challenge from numerical evidence of the presence of recurring transient 
spikes, which are non-local structures, in the approach to singularities [SI [7]. In- 
sufficiently resolved spiky structures in the singular regime have been observed 
in numerical simulations of C/(l) symmetric spacetimes [24], [171 Chapter 7] and 
generic spacetimes [T5] . 

In this paper I present the first explicit spike solutions in the class of Gowdy 
spacetimes. I then conjecture that spikes are an integral part of the dynamics 
on approach to generic singularities, and should constitute the non-local part of 
the generalized Mixmaster attractor. 
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2 Gowdy spacetimes 



Gowdy spacetimes refer to solutions of the vacuum Einstein equations, with 
metric of the form 

ds^ = _e(^-3-)/2dr2 + e(^+^)/2dx2 ^ ^^-^dy + Q dzf + c-^-^dz^. (1) 

The vacuum Einstein equations imply that P{t,x) and Q{t,x) are determined 
up to a constant by 

Prr = e'''Ql + e-'^P,,-e^''Ql) (2) 

Qrr = -2PrQr+C-^''iQ.. + 2P^Q^), (3) 

while A(r, a;) decouples from the above equations and is determined up to a 
constant by 

A. = -{P^ + e'^Ql) - e-'^P'^ + e^^Ql) (4) 
K = -2{PrP.,+e^PQ,Q,). (5) 

The decoupling of A simplifies the analysis of Gowdy models, as focus can be 
narrowed down to P,Q. The time variable t tends to infinity as the singularity 
is approached. The so-called G2 spacetimes are those which admits two com- 
muting Killing vector fields acting in the (y, z) planes (or cylinders or tori), and 
Gowdy spacetimes are a special case, in which the action of the G2 group is 
orthogonally transitive. 

The dynamics of Gowdy spacetimes is of interest as part of the larger set of 
oscillatory behaviour, because it describes the dynamics during a Kasner era (see 
e.g. [m Section 7]). Although the Kasner era in Gowdy spacetimes terminates 
at Kasner solutions, in the more general G2 models a Kasner era is followed by 
another Kasner era, without termination, hence giving the oscillatory Mixmaster 
dynamics. 



3 Solution-generating transformation 

The solution-generating transformation in interest is essentially the one used in 
[4]. It is composed of two transformations. The first is modified inversion in the 
(P,Q) hyperbolic plane [4],[T1 eq (3.11)]: 

where I have put a minus sign in Q to make it a reflected inversion, so that it has 
a better interpretation as a frame rotation below. A is unaffected. The second 
is the solution-generation transformation called the Gowdy-to-Ernst transfor- 
mation [1]: 

P = -P + T, = -g2(^-^)q„ Qx = -c'''Qr. (7) 
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The obstacle in generating explicit solutions is the integration to obtain Q. 
Composing the two transformations gives a powerful solution-generating trans- 
formation that can be iterated to generate a family of new solutions. 

The transformation and the solutions can be presented more elegantly in 
the orthonormal frame formulation with the so-called /3-normalized variables as 
presented in [TBI Section 4.4.1], in the so-called timehke area gauge with 

t = T, A/o-i = -2, Ei^^2e-\ (8) 

In the orthonormal frame formulation, an orthonormal frame is used, and vari- 
ables are the frame components and (essentially) their first derivatives, divided 
by /3. The {y, z) area expansion rate (3 is related to A by T^, App A. 3] 

/3 = _ie-(^-3-)/4. (9) 

Negative value for /3 describes contraction as r increases. 

For Gowdy models, the key /3-normalizcd variables arc (S_, A^x , Sx , iV-), 
which are decompositions of orthonormal frame components of the 3-by-3 Sa/j 
and Naf3 matrices: 



^cp^ \ S+ + V3^- \/3Sx (10) 




A^„^ = I 2V3N^ VSNy, I , (11) 
where E+ is given by 

E+ = i(l-I]2_ -Ar2 -7V2). (12) 

In order to preserve the form of the spatial orthonormal frame itself is 

rotating around the x-axis at the rate of i? = — VSSx [TB] Section 3.1]. Under 
this condition, (S_, A^x , Sx , A^-) are related to the derivatives of P, Q by 

S- = -^, A^x=-^^, Ex=-^^, N^^^^-^. 
V3 V3 V3 V3 

(13) 

The reflected inversion transformation has a simple interpretation as the 
rotation of the spatial orthonormal frame to another rotation rate that also 
preserves the form of (fTT|) . The angle of rotation (p is determined by 



(Qe^)2-1 . ^, 2Qc^ 



(E_,Ex) rotates as follows: 



E_\ _ / cos2(/) sin2(/«\ /S_ 
Sx/ ^ V-sin20 cos20) IEx 



(15) 
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and similarly for {N^ , ) . 

The Gowdy-to-Ernst transformation is much simpler in terms of , A^x , , -/V_ 

(S_,7Vx,Sx,iV_) - - ^,-iVx,iV_,Sx). (16) 



4 The explicit solutions 

I now present the Kasner solutions, the rotated Kasner solutions, the Taub 
vacuum solutions, the rotated taub solutions, and the spike solutions. I shall 
discuss their dynamics in the next section. 



4.1 The Kasner solutions 

We shall use the Kasner solutions (vacuum Bianchi type I solutions) as the seed 
solution. The metric components (P, Q, A) are given by 

p = iuT + Po, Q = Qq, A = -wV + Ao, (17) 

with w serving to parametrize the Kasner solutions, w is related to the Khalatnikov- 
Lifshitz parameter u hy w = 2u + 1. Pq, Qq and Aq are arbitrary constants and 
do not parametrize the Kasner solutions. 

In terms of /3-normalized variables, the Kasner solutions are given by 

11) 

(S_,iVx,Sx,^-)- (-^,0,0,0). (18) 

The flat Kasner solutions are those with w = ±l,oo, and the other plane- 
symmetric Kasner solutions are those with w = 0, ±3. 



4.2 The rotated Kasner solutions 

It is essential to choose Qo 7^ so that the rotation transformation is nontrivial. 
Applying the rotation transformation ([6]) to the Kasner solutions yields the 
Kasner solutions viewed in a rotating frame (see also [U eqs (3.12)-(3.13)], O 
eq (7)]). The metric components (P, Q, A) are given by 



e 



-p 



-{WT + Po) 



Q2+e-2{wr+Po) (20) 

A = -wV + Ao. (21) 

It turns out to be convenient to write P and Q in the form 

p sech(ii;T + Pn + In Qo) ^ 1 r , / „ , ^ m ^„„x 

e-^ = ^ — ^ Q = - — [l + tanh(wT + Po+lnOo)]. (22) 



(19) 
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For simplicity wc will now choose Pq + InQo = 0. The /^-normalized variables 
are then given by 

( w w \ 

(S_,iVx,Sx,Af-) = f --^tanh(wT),0, -^sech(wr),Oj . (23) 

4.3 The Taub vacuum solutions 

Applying the composed transformation ©-Q on the Kasner solutions with 
w; ^ (For the w = Kasner, the transformation gives the w = \ Kasner.) 
yields the Taub vacuum solutions (vacuum Bianchi type II solutions), whose 
(P, Q, A) are given by 

P = r + ln(sech(wr))-ln(2Qo) (24) 
Q = 2Qowx + gi (25) 

A = -2 ln(sech(u;r)) - (v? + l)r + Ai. (26) 

Qi and Ai are arbitrary constants and do not parametrize the Taub solutions. 
The /3-normalized variables are given by 

(S_,iVx,Ex,^-) = (^-^tanh(z«T)--^,0,0,-^sech(«;r)j . (27) 

4.4 The rotated Taub solutions (the false spike solutions) 

Note that the frame of the Taub solutions is not rotating. Applying the rotation 
transformation ((6|) to the Taub solutions yields the rotated Taub solutions: 

P = -r - ln(sech(wT)) - ln[(u;e^sech(wr)x)2 + 1] + ln(2Qo) (28) 

Q ^ 1 wa;(e^sech(wr))^ 

^ 2Qo {we^sech{wT)xy + 1 ^ ^ 

A = -2 ln(sech(u;T)) - {w^ + 1)t + Ai, (30) 

where for simplicity we have chosen Qi = 0, so that the term x — Qi/ {2wQo) 
simplifies to x. The /3-normalized variables are given by 

(E_,iVx,Sx, A^-) = (cS-Taub, -SA^- Taub, -sS_Taub,cA_ Taub) , (31) 

where 

P-1 2/ 



f = (Qe )Taub = u;e'^sech(wT)x. (32) 



The rotation transformation makes the frame rotates in a spiky manner near x = 
0. The rotated Taub solutions are also referred to as the false spike solutions. 
The shape of the spike will be described in the next subsection. 



5 



4.5 The spike solutions 

Applying the composed transformation ©-Q on the Taub vacuum solutions 
yields a new solution, whose (P, Q, A) are given by 

P = 2t + ln(sech(wr)) - ln[(we^sech(wT).T)2 + 1] - ln(2Qo) (33) 

Q = -4Qow[e"^^ + 2(wtanh(wr) - l)x^] + Q2 (34) 

A = -4 ln(sech(u;r)) + 2 ln[(we^sech(wr)a;)^ + 1] ~ (w^ + 4)r + A2. (35) 

The /^-normalized variables are given by 

(S_,Afx,Sx,iV_) = ^-cS_Taub- ;^,sA''_ Taub, cA''- Taub, -sS.Taub^ • (36) 

The spatial dependence of (E„ , A^x , Sx , A^-) lies in c and s, which depends on 
X in a spiky way. c has a spike of the form 

which has a single peak at a; = and two zeros at a; = ±l/fc. s has a spike of 
the form 

2kx 

{kxY + r ^^^> 

which has two peaks with opposite signs at x — ±l/k and a zero at a; = 0. From 
((32|) . we have k — we'^sech{wT), so the location of the peaks varies with time. 
The limit of fc as r — > oo bifurcates at \w\ = 1, as can be seen clearly when k is 
written in the form 

fc = u'e^-l'"^l(l +tanh|wr|). (39) 
For \w\ > 1, fc is largest (and the spike is narrowest) when 

r = — arctanh— , (40) 
w w 

and k tends to zero as t ~f oo. For < jwl < 1, fc tends to infinity (and the 
spike increasingly narrow) as r — )■ oo. 

4.6 The radius of the spike 

Therefore, we shall define the radius of the spike to be one half the x-coordinate 
distance between the two peaks in s, which works out to be 

radius of spike — — c^^ coshwr. (41) 

w 

We must measure the radius of the spike against another important scale - the 
particle horizon (when approaching the initial singularity) or the event horizon 
(when approaching the final singularity), whose x-coordinate radius is given by 

/oo 
IMolEi^dr = e~'^, (42) 
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Figure 1: The spiky structures of the s and c functions in ((32)) with w — 1 
against r and X = e^x. 



as follows from ([8|). Thus, the radius of the spike as a multiple of the horizon is 



This measure of the radius is more meaningful than (|4T|) . Using (|43|) . the spike 
is narrowest relative to the horizon at r = 0, which does not coincide with (|40|) . 

Notice that for < |w| < 1, the size of the spike is always super-horizon^ 
despite appearing to be increasingly narrow. On the other hand, for \'w\ > 
1, there is a brief period during which the spike spike becomes a sub-horizon 
structure. Physically, during this period, an observer along the spike worldline 
a; = would see substantial inhomogeneity. Mathematically, the second-order 
spatial derivative terms P^x and Q^x in the Einstein equations are not negligible 
during this period. 

To see the sub-horizon structures, it is useful to plot variables against mul- 
tiples of the radius of horizon, i.e. against the variable 



in which the horizon of the observer along a; = is represented by X ~ ±1. 
Figure [T] plots the spiky structures of the s and c functions in ([32]) against r 
and X. 

5 Visualizing the dynamics of the spike solu- 
tions 

To visualize the dynamics of the spike solutions, we will use both (3- and Hubble- 
normalized variables. For the spatially homogeneous background dynamics, it 
is best to use the Hubble- normalized variables (see [20l Chapter 6]), which are 



radius of spike — — cosh wt x radius of horizon. 

w 



(43) 



X = e'x. 



(44) 
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w=-3 




w=3 



w=1 



Figure 2: The Kasner circle and the parameter w. 



related to the /S-normalized ones via 



(S+,S_,Ex,iV_,iVx)^ 



1 - 



-(S+,S_,Sx,iV_,iVx), (45) 



and satisfy 



^H2 



^H2 



^H2 



H2 



N"^ = 1. 



(46) 



In the state space of Hubble-normalized variables, the Kasner solutions appear 
as equilibrium points on a unit circle in the (S]_|_,E_)^ plane (see Figure [5]). 
The three different representations of the flat Kasner solution are marked by 
a circle, and the three representations of the other axis-symmetric Kasner so- 
lution are marked by a triangle. All other Kasner solutions have six different 
representations, each located equidistant from a flat Kasner point. 

Each ui-rotated Kasner solution appears as a semi-circular orbit on the 
sphere T,^"^ + T,^"^ + — 1 (Figure [3]), connecting two representations of 
the same Kasner solution. As r goes from —oo to oo, the w-rotated Kasner 
solution tends from the — |?i;|-Kasner solution to the |?ii|-Kasner solution (also 
compare the limits of ([23| with p^ ). The sign of w determines the sign of E^. 
The set of all such orbits (all values of w 7^ 0) is referred to as a frame transition 
set [H]. 

Each w-Taub solution appears as a semi-circular orbit on the sphere Y,^'^ + 
Y^H2 _^ ]\jH2 _ 2 (Figure lU, connecting two Kasner points. As r goes from 
—00 to 00, the w-Taub solution tends from the (\w\ + 1)-Kasner solution to the 
(1 — |w|)-Kasner solution. The sign of w determines the sign of . The set of 
all such orbits (all values of w 7^ 0) is referred to as a curvature transition set. 

When projected on the (I]+,I]_)^ plane, the frame transition sets form 
parallel straight lines, while curvature transition sets form straight lines em- 
anating from one corner of a triangle superscribing the Kasner circle (Figure 
[5|). The two transition sets combine to describe the dynamics during a Kasner 
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Figure 3: The rotation orbits (frame transition set) on the sphere S;^^ + S;?^ + 
= 1. Arrows indicate the direction of increasing r, towards the singularity. 




Figure 4: The Taub orbits (curvature transition set) on the sphere 



H2 



^H2 



era, which consists of long Kasner epochs (described by the Kasner equilibrium 
points), punctuated by brief periods of transitions. During frame transitions, 
the rotating frame interrupts a Kasner epoch with a brief period of frame rota- 
tion, thus fictitiously splits the same Kasner epoch into two epochs - one with 
negative w and one with positive w. During curvature transitions, the spatial 
curvature becomes significant. Successive pairs of curvature and frame transi- 
tions reduce the value of the Kasner epoch by 2, and terminate at the final 
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Figure 5: The rotation and Taub orbits projected on the plane. 



Kasner epoch with < w < 1. For examples, 

u; = 4.2 -2.2 2.2 '-^ -0.2 0.2, (47) 
^ 5.2 -3.2 3.2 -1.2 1.2 "-^ 0.8. (48) 

See [T^ for a detailed discussion of these solutions. 

Applying the rotation transformation to the Taub solutions yields the false 
spike solutions, which are simply the Taub solutions presented in a spatial frame 
that rotates in a spiky way. The Gowdy-to-Ernst transformation then maps the 
false spike solutions to the real ones. The dynamics of the false spikes is similar 
to that of the real spikes (related through the simple map (fT6|) ). so I shall focus 
on the real ones. 

For the spike solution with value w, as r goes from — oo to oo, all orbits 
begin at the same Kasner point (with value \w\ + 2). For \w\ > 1, the orbits also 
end at a common Kasner point (with value 2— The smooth limit describes 
the transient spikes. For < |w| < 1, however, the orbit along x — ends at 
the Kasner point with value 2 — while orbits along x ^ ends at the Kasner 
point with value \w\. This discontinuous limit describes the permanent spikes. I 
shall discuss this in more detail in the next subsection. The = 1 bifurcation 
of the limit is due to the limit of the factor sech('u;r)e'^ as r ^ cx) (see eq (P^ ). 

Along worldlines far away from the spike worldline, the orbits approximate 
the Taub and rotation orbits. Along the spike worldline, the orbit (called the 
spike orbit) lies on the sphere S;^^ + S;?^ + = 1, and forms a straight 
line when projected on the (S+,S_)^ plane. Along worldlines near the spike 
worldlines, the orbits interpolate between the two extremes. The set of all 
orbits (all values of w ^ 0) is called the spike transition set. Figure [S] shows 
four members of the spike transition set. 
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w=0.5 




Figure 6: Orbits of the spike solutions projected on the plane. Orbits 

are coloured red (thick arrowed line) along the spike worldline a; = 0, blue (thin 
arrowed line) along x = 1000, and magenta (thin line without arrows) along 
small values of x. w = 0.5, 1, 1.5, 3.5 respectively. 



The inhomogeneous spike dynamics are best visualized in the /3-normalized 
N- and Ex, which are active variables in the Taub and rotation orbits respec- 
tively, and considered the active players in spike dynamics. Figure [7] shows the 
spikes in /3-normalized A^- and Ex- 
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5.1 Transient and permanent spikes 

We shall discuss the dynamics of transient spikes first. For a spike solution with 
jwl > 1 one gets a transient spike that smooths out as t ^ cx3. In this case 
there is an instant in time (which coincides with (|40p ) 

Tflip = — arctanh— , (49) 

WW 

when the foUowing variables become homogeneous: 

1 3 — w^ 

^-(Tfiip,a:) = 0, E_(rflip,a;) = i;+(rfiip, x) = — - — . (50) 

The spike in N- flips sign at r = Tflip. All orbits also project on the same 
point on the (I]+,I]_)^ plane at this instant in time. See Figures [6] and [71 
Dynamically, faraway observers undergo three transitions - curvature, frame 
and curvature: 

I I „ curv I I frame , , curv „ , , 

|u;| + 2 — > ~\w\ — > \w\ — >2-\w\, (51) 

where the values |w|+2, —\w\, \w\ and 2 — are those of the consecutive 
Kasner epochs in the spike solution. While faraway observers undergo the first 
curvature transition, the spike observer, because (r, 0) = 0, is unable to do so, 
and observes the formation of a spike, which starts as a super-horizon structure. 
While faraway observers undergo the frame transition, the spike has narrowed 
to become a sub-horizon structure, reaches a minimum width at r = 0, and 
then starts to widen, and flips sign at r = Tflip. The spike observer transitions 
directly to the fourth Kasner epoch: 

\w\+2 — >2-\w\. (52) 

While faraway observers undergo the second curvature transition, the spike 
continues to widen and eventually smooths out. 

We now discuss the formation of permanent spikes. Recall that Gowdy 
models can undergo only one Kasner era. The w-valuc of the final Kasner 
epoch of the era satisfies < w < 1. There are two different penultimate 
Kasner epochs - one with 1 < w < 2 and the other with —l<w<0. The 
former transitions to the final Kasner epoch through the curvature transition, 
and the latter through the frame transition. See examples (|47| - (|48p . 

Permanent spikes form in the spike solutions with < |w| < 1. Dynamically, 
faraway observers undergo two transitions - curvature and frame, and terminate 
at the final Kasner epoch: 

\w\+2^-^ -\w\'°'-^°''\w\. (53) 

The spike observer on the other hand transitions to the penultimate Kasner 
epoch with 1 < w < 2 as its final Kasner epoch: 

|w|+2 — y2~\w\. (54) 
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x10^ xlO* 




Figure 8: The Weyl scalars for the spike solution with w = 2.5. 



The resulting discontinuous pointwise (fixed x) limit is the permanent spike. 
One at first gets a transient spike that forms during the antepenultimate Kasner 
epoch with 2 < w < 3, that later becomes a permanent spike after the transition 
((M)) . See Figure [6] with w = 0.5. 

Unfortunately, the spike solutions do not include permanent spikes that form 
directly during the penultimate Kasner epoch with 1 < w < 2. Spikes that form 
during the other penultimate Kasner epoch with — 1 < w < are false spikes. 
Spikes do not form during the final Kasner epoch (with < w < 1). 

5.2 The Weyl invariants 

For the spike solutions, the explicit expressions of the Weyl invariants (formulae 
given in Appendix), are too complicated to be useful. Nonetheless we can 
compute them numerically to tell apart real and false spike solutions. The Weyl 
scalars for a spike solution are shown in Figure [51 Observe the transient, sub- 
horizon structures in the Weyl scalars. The Weyl scalars blow up as t — > oo, 
indicating the approach to a curvature singularity. The Weyl scalars for the 
Kasner solutions, the rotated Kasner solutions, the Taub solutions and the false 
spike solutions are spatially homogeneous. 
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w=0.75 



w=1.25 




Figure 9: Two families of orbits with w = 0.75 and 1.25. The Kasner seed is 
indicated by a black *, the rotation orbit in light blue (thick dashed line), the 
Taub orbit in dark blue (thick solid line), orbits of the false spike solution in 
magenta (thin dashed line), and orbits of the spike solution in red (thin solid 
line). 



5.3 Spike as the non-local part of the generalized Mix- 
master attractor 

The spike solutions closely resemble those produced numerically (see e.g. [5l[2])- 
This raises the question whether the spike solutions arc an integral part of the 
generalized Mixmaster attractor. 

Applying the rotation transformation and the Gowdy-to-Ernst transforma- 
tion to the Kasner solutions successively generates solutions describing the frame 
transition, the curvature transition, the false spike transition, and the (true) 
spike transition. See Figure [5] for two families of orbits generated. The Kasner 
circle, the frame transition set and the curvature transition set arc known to 
be part of the generalized Mixmaster attractor (see e.g. 14, Section IV], J3]). 
Thus it is natural to expect that the transformations preserve this property, 
and extend it to the false and true spike solutions, provided that spikes do not 
disappear. This is supported by numerical evidence, which indicates that spikes 
generally recur on the same spot (where N- = 0), despite numerical evidence 
that pairs of zeros of can be created or annihilated. See [151 concluding 
remarks] for a discussion of interaction between spikes. 

I therefore conjecture that the spike transition set (as well as the false spike 
transition set) constitutes the non-local part of the generalized Mixmaster at- 
tractor. The local part of the BKL conjecture, which states that the dynamics 
becomes asymptotically local as the singularity is approached, should be modi- 
fied to allow for spikes, whose dynamics is non-local, but only occurs within the 
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horizon of isolated worldlines. 



5.4 The simplest gauge for the spike solutions 

A natural question arises regarding the merit of using a rotating frame, which 
produces false spikes. It turns out that the spike solutions are not elementary in 
a non- rotating frame (because the time integral of Ex in ([36]) is not elementary). 
Furthermore, the net rotation in the spike solutions is non-zero and depends on 
X. As a result, the spike orbit does not appear as a straight line when projected 
on any plane, and all other orbits end on Kasner points lying on the Ex)^ 
plane with fixed E;^. Therefore a simple description of the spike transition set 
in a rotating frame is preferable to a complicated one in a non-rotation frame. 

We now consider other temporal gauges. The flipping of the spike occurs 
instantaneously in the timelike area gauge [TB], but not in other gauges. This 
singles out the timelike area gauge as the simplest temporal gauge to present 
and study the spike solutions. 

5.5 Higher-order spike solutions 

Applying the transformation on the spike solutions yields another new solution, 
whose (P, Q, A) are given by 

P = 3r + ln(sech(wT)) - \n[{we' sec\i{wT)xf + 1] - ln[/| + 1] - ln(2Qo) (55) 

Q ^ -iQoe^'x 3[(u;e^sech(wT)a;)2 + 1][(/| - l)w + 2/2(wtanh( wt) — 2) coshuirl 

-|-4we2^a;2(w2 + 2-3wtanh(wT)) e"^^ + Q3 (56) 

A = -61n(sech(wr)) 41n[(we^sech(wr)a;)2 + 1] + 21n[/| + 1] - (w^ + 9)r + A3, 

(57) 

where /2 is the factor Qe^ of the spike solutions. These second-order spike 
solutions have multiple zeros of N_ within the horizon of a; = 0. The solutions 
are parameterized by w and Q2/Qo- 

The variables (E_, TVx , Ex , A^-) of the third-order spike solutions can be 
obtained. A_ has even more zeros within the horizon of a; = 0, These solutions 
are parametrized by w, Q2/Q0 and Qa/Qo- 

These second and third-order spike solutions have not been commonly seen 
in numerical simulations because initial data with multiple zeros of A_ within 
one horizon are not commonly chosen in numerical simulations. 

It is of interest to study the movement of these zeros over more than one 
Kasner era, in the context of G2 models. If multiple zeros stay within the 
same horizon, then high-order spike solutions are conjecture to be part of the 
generalized Mixmaster attractor. Conversely, if they eventually move away from 
each other's horizon, then all high-order spike solutions gradually splinter into 
first-order spikes, and are not part of the attractor. I shall leave the detailed 
analysis of higher-order spikes for future works, and note its potential relation 
with the high- velocity spike approximation in ^ Section 6] . 
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5.6 Other seed solutions 



Seed solutions other than the Kasner and Taub solutions can be used. The 
following are examples of explicit solutions in the class of Gowdy spacetimes. 

• Solutions with Q = const and separable P — F{t)G{x) and their linear 
combinations [HI Section 3]. 

• Solutions (apparently new) with Q = const and additively separable P: 

P ^ ci{2x'^ + e-'^^) + C2T + C3X + a. (58) 

• The Wainwright-Marshman solution [551 Case I with ni = —j^], [71 Sec- 
tion 6.2]. 

But because the Kasner and Taub solutions are part of the generalized Mixmas- 
ter attractor while these other solutions are not, any spike solutions generated 
from these solutions are not expected to be part of the attractor. 

6 Conclusion 

Gowdy spacetimes have been appreciated as a simple class of inhomogeneous 
spacetimes for both manageable mathematical and numerical analyses. Even 
so, the discovery of the spike solutions is a pleasant surprise. 

The main results of this paper are the derivation of the explicit spike solu- 
tions using a solution-generating transformation, and the conjecture that the 
spike solutions are part of the generalized Mixmaster attractor. Higher-order 
spikes remain to be analyzed, and may also be part of the attractor. 

A similar solution-generating method was used in [23] to generate exact 
inhomogeneous U{1) symmetric solutions with Mixmaster dynamics. Although 
these solutions contain G2 solutions as a special case, comparison with the 
spike solutions is difficult due to the difference in the temporal gauges used. 
Nonetheless, it is important to examine if these solutions also contain spikes. 

It is reassuring to see the spike solutions very closely resemble those seen 
numerically. The spike solutions will be immensely valuable to the study of 
spikes in more general classes of spacetimes. I expect the spike solutions in this 
paper to be only one of several possible kinds of spikes in generic spacetimes. 
Insufficiently resolved spiky structures in the singular regime have been observed 
in numerical simulations of C/(l) symmetric spacetimes [53],[T71 Chapter 7] and 
generic spacetimes [18]. Therefore, the ultimate goal is to understand these 
spikes, and to determine whether they are part of the generalized Mixmaster 
attractor in generic spacetimes. Numerical study of spikes in G2 cosmologies in 
collaboration with David Garfinkle, Frans Pretorius and Lars Andersson is in 
progress. 
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A The Weyl scalar invariants 

The orthonormal frame components Catcd of the Weyl tensor can be conve- 
niently expressed in terms of the electric and magnetic components Eafs and 

CaOpn — Eal3, Ca/SjS — —^^a^^S^l^v, Ca/SjO — a/sHfii, (59) 

which are then normalized by 3/3^: 

and further decomposed the same way as the shear matrix. The components 
are given by 



are computed as follows 









i(l-3E+)E_ + . 


fx - 


i(l-3E+)Ex +■ 


n+ = 


-N^H^ -iVxSx 


H- = 


-S+Ar_ - fiV+E. 




-E+iVx - f iV+S 


V3^_ 


, r = -3(iVxE_ - 



l{Nl+Nl] 



(61) 



3-+— ■ W. - '^)^x (62) 



dx 
d_ 
dx 



N_ (63) 
(64) 

i(2e-A_,)Sx (65) 
i(2e-A_,)s_ (66) 



CabcdC^'^" = KEcpE^^ - H^pH^P) (67) 

CabcdC^''"' = l(!>Ec.pH^^ (68) 

Cab^^C.d'fCef'''' = -16(£;/£;/3^i?7 - SEjHp^H^"^) (69) 

Cab'''C,f^*Cef'"' = leiHjHp^H^'' ~ SEjEp^H^''), (70) 

where *Cabcd = ^Vab'^'^ Cefcd, f-nd Tj"'^'''^ is the totally antisymmetric permutation 
tensor, with 7^0123 _ ^ 
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